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COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


THIS SOFTWARE s FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ge VE COPYRIGHT NOTICE. THIS SOFTWARE OR ANY OTHER 
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FACILITY: MATH LIBRARY 
“ABSTRACT: 


MTHSHLOG and MTH$HLOG10 are functions which return the H floating point 
natural or common logarithm of their H 6 Ra and peice argument. The call 


00 is standard call-by-reference. BL... HLOG_R8 an Senger 0_R8 are spec tet 
0000 routines which are the same as MTHSHLOG “and MTHSHLOGIO éxcept 
0000 faster non-standard JSB call 13 used with the sraaaal in RO through R3 
4 and no registers are saved. 
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AUTHOR: 
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Edit History for Version 1 of MTHS$HLOG 


1-004 - Adapted from MTHS$GLOG version 1-001. JAW 24-Sep-1979 
- eae geete rene after POLYH for greater accuracy. JAW 
6- n-1 
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; - 01 - Add MTHSHLOG2. RNH 08-Aug-1981 
; Correct %. lo Is, al JSB entries. Use G* addressing for 
externals. SBL apt 
erenged MTHS$AB ALOG © MTH$$AB_ALOG_V. RNH 29° lade 
nated symbolic short Literals. “RNH 15- 1 
Changed H_FHI to the global symbol MTHS$$AB_H_ FHI. PDG 3-Nov-81 
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-SBTTL DECLARATIONS 3; Declarative Part of Module 


$ INCLUDE FILES: MTHJACKET.MAR 

>: EXTERNAL SYMBOLS: 
SEXTRN lhe oe 3; Error code 
SEXTRN S1G + Math $i nal routine 
eEXTRN NTHSSAB” CRLOG. Vv : Table o byte offsets 


EQUATED SYMBOLS: 


ACMASK = “M<IV, R2, R3, R4, RS, R6 R7, R8&> 
 heatster | save mask and IV enable 


MACROS: none 
PSECT DECLARATIONS: 


-PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 
$ Program section for math routines 
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put te bECLARAT ION poing Natural ang Com ) "Modul "BrSEb=1 $84 i 38: 33 UATARTL. SREIATHRLOG. MAR; 1 _ 5) 1 
DDA8 SFFF -OCTA “XODA8 9 SFFF 3C2 = .4000000000000000000000000034 
SBEs GSuG6SGESSHGS GERSRS HESS S86 REA SRR start EDS SSHISHISHI SSNS 
0000 3 ze | LOGLEN2 = .-LOGTAB2/ 
re 26 H_LN_2_HI1: 
00006730 93C7F357 AS9E2FEF 62E44000 22 28 ey Taha “x0000673093C7F 357A39E2F EF 62E44000 ; Hi 98 bits of Ln2 
069E16C5 4C€5B89339 79A157A0 F9O7B3F9A rf rs = snase OCTA “X069E16€54C58933979A15 het 4) rhe ; Low bits of Ln2 
6508152 BCB7SE FF a670 6g ~ “LONG SXBCB73EEE. *x66508152 + 3e134$5448190325182765112891891660508 
tt 64 H_INV_ Une, CONS: ‘ 2 + convert from natural log to log base 
B82F7652 71544001 0680 65 UAD XB82F 765271544001 sé 
D25AFDAO 7D0FE177 Beas 66 “SUAD “XD25AFDAO7D0FE177 
0690 68 
0690 269 
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ndard H-Floating LOG SEP=19 MTHRTL.SRCJMTHHLOG.MAR; 1 
-SBTTL MTHSHLOG = Standard H-Floating LOG 


t 
0 


eo 


fthsa 
: FUNCTIONAL DESCRIPTION: 

; HLOG =H floating point LOG function 

: HLOG(X) is computed using the following approximation technique: 
If X =< 0, error. Otherwise 

Let X = f * (2**#n), where 1/2 <= f < 1 


If n is yrgeter than or equai to 1 than 
et set N=n- 1 and F = 2ef. 
ge 


set N =n and F = f. 
Then Ln (x) = Neln2 + Ln(F) 
If iF = 13 < 2e*-5 then 

in 


(F) = W + WeP(W), where W = F - 1 and P 
is a polynomial of degree 18. 
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0690 Else 
0690 Ln(F) = Ln(FHI) + Z*Q(2*Z), where FHI is ob- 
0690 tained by table look-up, Q is a Pel ynentes of 
+24 degree 10 and Z = (F = FHI)/(F + FHI) 
0690 NOTE: The quantities Ln(FHI) and ln2 are used in the above 
0690 equations in two parts - a high part (containing the 
0690 h . order bits) and a low part (containing the low 
0690 order bits. In the code the high and low parts of the 
0690 constants are indicated by a _HI and _LO suffix respec- 
0690 tively. The values were chosen such that N*LN_2_HI1 + 
0690 LN_FHI_HI is exactly representable. 
0690 ; 
0690 ++ 
0690 : 
0690 ; CALLING SEQUENCE: 
0690 10 ; 
0690 1 hlog.wh.v = MTHSHLOG(x.rh.r) 
0690 1 -or- 
8098 } CALL MTHSHLOG(hLog.wh.r, x.rh.r) 
$690 15 ; INPUT PARAMETERS: 
0690 16 at 
00000004 0690 17 LONG = 4 ; Define Longuord multiplier 
00000008 24 3 x = 2 * LONG : Contents of x is the argument 
£30 9 ; IMPLICIT INPUTS: none 
o88 ¢ 3; OUTPUT PARAMETERS: 
00000004 $30 : y hlog = 1 * LONG ; Contents of hlog is the result 
4 § 3 VALUE: H floating logarithm of the argument 


— 
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$3 8 : IMPLICIT OUTPUTS: none 

69 0 : COMPLETION CODES: none 

69 : : SIDE EFFECTS: 
0690 4; Signals: MTH$_LOGZERNEG if !X! =< 0.0 with reserved operons in RO/R3 (copied 
690 5 ; to the signal mechanism vector CHFS$L_MCH_RO/Ri by LIBSSIGNAL). Associated 
690 6 ; message is: ‘‘LOGARITHM OF ZERO OR NEGATIVE VALUE’’. Result is reserved operand 
$38 38 ; “0.0 unless a user supplied (or any) error handler changes CHFS$L_MCH_RO/R1. 
+24 § : NOTE: This procedure disables floating point underflow, enables integer over- 
69 40 ; flow, causes no floating overflow or other arithmetic traps, and preserves 
+24 41 : enables across the call. 

069 re : 

0690 43 ;--- 
0690345 

41FC 0690 346 -ENTRY MTHS$HLOG, ACMASK ; Standard call-by-reference entry 
Be3s 34 ; Disable DV (and FU), enable IV 
be8 348 MTHSFLAG_JACKET ; Flag that this is a jacket procedure 
6D Q0000000'GF 9E $693 MOVAB G*“MTHSSJACKET_HND, (FP) 
0699 ; set handler address to jacket 
994 ; handler 
0699 349 : in case of an error in special JSB 
0699 350 3; routine 
50 08 BC 70FD 0699 32) MOVH ax(AP), RO ; RO/R3 = arg 
44 069E 26 BSBB MTHSHLOG_R8 ; Call special HLOG routine 
04 BC 50 7DFD O6A0 35 MOVO RO, @hlog(AP) ; Store result in first argument 
04 se $2¢ RET ; Return to caller 


6D 


00000004 
00000008 


00000004 


41FC 


QO000000'GF YE 


50 
04 BC 


08 BC 70FD 
if 10 

50 7DFD 

04 
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-SBTTL MTHSHLOGIO = Standard H Floating Common Logarithm 


FUNCTIONAL DESCRIPTION: 
HLOG1O - H floating point LO0G19 function 
HLOG10(X) is computed as HLOGIO(E) * HLOG(X). 
See description of MTHSHLOG 
CALLING SEQUENCE: 
hlog_base_10.wh.v = MTHSHLOG10(x.rh.r) 
CALL MTHSHLOG10(hLog_base_10.wh.r, x.rh.r) 
INPUT PARAMETERS: 


; Define me multiplier 
; Contents of x is the argument 


LONG = 4 
x = 2 * LONG 

OUTPUT PARAMETERS: 
hlog_base_10 = 1 * LONG 


SIDE EFFECTS: See description of MTHSHLOG 


-ENTRY MTHS$HLOG10O, ACMASK § siendere call-by-reference entry 
; Disable DV (and FU), enable IV 
MTHSFLAG_JACKET ; Flag that this is a jacket procedure 


MOVAB G*MTHSSJACKET_HND, (FP) . 
; set handler address to jacket 
; handler 


in case of an error in special JSB 


RO/R3 = arg ; 

Call special HLOG10 routine 
Store result in first argument 
Return to caller 


ax(AP), RO 
BSBB MTH$HLOG10_R8 
aoe RO, @hlog_Base_10(AP) 
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ets : 9 -SBTTL MTHSHLOG2 = Standard H Floating Common logarithm | 
6BC 402 ;++ 
re$3 : 3; FUNCTIONAL DESCRIPTION: | 
3 2 5 Z HLOG2 =H floating point LOG2 function 
b6B¢ é : ; HLOG2(X) is computed as HLOG2(E) * HLOG(X). 
peat 609 ; See description of MTHS$HLOG 
b6B¢ 41} ; CALLING SEQUENCE: 
Rees “i$ : hlog_base_2.wh.v = MTHS$HLOG2(x.rh.r) 
: -or- 
poet 212 5 CALL MTH$HLOG2(hlog_base_2.wh.r, x.rh.r) 
d6B¢ 617 + INPUT PARAMETERS: 
00000004 068C 419 © LONG = 4 ; Define longword multiplier 
00000008 ayes r Y x = 2 * LONG ; Contents oF x is the praument 
06BC 422 : | 
pepe ? ; OUTPUT PARAMETERS: 
00000004 06B¢ 4 Bs hlog_base_2 = 1 * LONG | 
06BC 4 5 : SIDE EFFECTS: See description of MTH$HLOG 
068C 428; 
068C 429 ;-- 
06B8C $39 
41FC $36 -ENTRY MTHS$HLOG2, ACMASK 3; Standard call-by-reference entry 
43 ; Disable DV (and FU), enable IV 
434 MTHSFLAG_JACKET ; Flag that this is a jacket procedure 


ooooo°coe 
SBEeeee 


6D 00000000'GF 9E MOVAB G*MTHSSJACKET_HND, (FP) 
; set handler address to jacket 
Bete ; handler 
06C5 435 ; in case of an error in special JSB 
06C5 $36 ; routine 
50 08 BC 70FD 06C5 43 MOVH ax(AP), RO ; RO/R3 = arg | 
O6CA 438 BSBB_ - MTHSHLOG_RB ; calculate natural log 
04 BC 50 BO AF 65FD 06C 4 MULH3 H_INV_LN2_CONS, RO, @hlog_base_2(AP) } 
06D 440 3; convert and store result in first | 
06D 441 3 argument 
04 0603 442 RET ; Return to caller | 
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on6 rt: -SBTTL MTHSHLOGHLOGIO_R8 = Special HLOG/HLOG10 routines 
ty 123 ; Special HLOG/HLOG10 - used by the standard routine, and directly. | 
604 468 ; CALLING SEQUENCE: 
604 449 ; save anything needed in RO:R8 
604 450; 4 Input in RO/R3 
0604 451: JSB MTHSHLOG10_ RB 8 /MTHSHLOG_ ke 
6D4 126 : return with result 
604 453 ; Note: This routine is aritten to avoid causing any ye ¥ overflows, 
604 454 ; floating overflows, or floating underflows or divide by 0 conditions, 
64 439 z whether enabled or not. | 
604 439 3 REGISTERS SED: 
604 458 ; RO/RS - H floating argument then result 
Bene rhe § eer - pareceeeteee results 
604 461; RB - Rothe wy into H_FHI table 
604 re | 
st ie | 
0604 465 MTH$HLOG10_R8:: : Spec tal HLOG10 routine | 
57 50 BO 0604 466 MOVW RO, R7 : = biased exponent 
08 15 06D7 467 BLEQ ERR : Error if <= 0 
0609 468 ; User PC on top of stack 
0609 469 : Note: ERROR rout ine B+ gue 9 on user 
0609 470 z PC being on top of 
res og) 3 used ne call to ATNSHLOR. R8 is not 
3 use 
OE 609 138 BSBB HLOG ne R8 ; Call_common Wi 06/HLOG10 routine | 
50 91 AF 64FD O6DB 474 MULH2 H_ LOG10_E, RO 3 RO/R3 = LOG10(e) * LOG(X) 
Reey $7? RSB 3; Return 
010A =«—i31 Bees pth ERR: BRW ERROR 
06E4 479 MTHSHLOG _RB:: 5 : ecial LOG routine | 
57 50 °9 06E4 480 Aovw RO, R7 3 = Biased exponent | 
F8 1 ed 4 rt 06 comin ae ERR 3 LOG (A) is not defined for x=<0 
57 4000 8F a 06E9 183 = $uBQ #*x4000, R7 3; R7 = Unbiased exponent 
6A 15 Q6EE 484 BLEQ NEG_EXP 3 Branch to processing for n=<0 
Bere a 
xponent is positive. =n-1an = 
60 2; E is positive. N 1 and F = 2f | 
6FO 489° 
57 B87 O6F 490 DECW R7 : - 1 
50 57 Ae F 491 SUBW R7, RO : RO/AS. = = 2f 
58 50 7 9C 6F £36 ROTL ae a0. R38 ; RB = he akg into MTH$$AB_ALOG table 
58 FFFFFFOO BF CA OGF9 49 BICL #- i a, 2 (o.exp bit and Ist'7 tract bits | 
56  00000000'GF oD 9 494 MOVAL G=ATWbSAB _ALOG_V, R6 =; - RO = Address of RIL vector ent | 
6 C ? 495 ADDL (R6) ; R6 = Address of MTH$S$AB_ALOG vod le 
99 70A 496 MOVB (Red ERBS R8 : R8 = offset into H_FHI tables 
1 7 43 BLSS LN_1_PLU Lug ; Branch to special processing 
71 498 : for F close to 1 
71 499 
71 500 
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; Floating Point Natural and Common 16-SEP-1984 01:36:4 AX/VMS Macro V04-00 Page 
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HH y ; Compute 2, Z2*#2, P(Z**2) and 2*P(Z2**#2) 
71 ~ 
7 57 GDFD 71 4 CVTWH R7, =(SP) ; Push N onto the stack 
58 FB8E6 CF4 FD 0714 05 MOVA MTHSSAB_H_FHICR8J], RB ; RB = Address of FHI 
54 63FD at $ SUBH (RB), RO, ~R4 + RG/R7 = F = FHI 
60FD 07 ADDH2 (RB)+, RO : RO/R3 = F + FHI 
ze 4 67FD 0724 8 DIVHS RO, R4, =(SP) : (SP) = 2 = (F = FHID/(F + FHI) 
5 E gf $3FD 729 «50 MULH3 (SP), (SP), RO + RO/R3 = Z*#2 
FE6B CF OA SFD (07 : 10 POLYH RO #LOGLENZ-1, LOGTAB2 : RO/R3 = P(Z*#2) 
0 BE 64FD , 3 i MULH2 (SP)+, R : RO/R3 = 2*P(Z*#2) 
739 13; 
, 3 1 ; Compute B = N*LN_2_LO + LN_FHI_LO + Z*P(7*2) 
54 FF21 CF 3 65FD 0739 1 MULH3 (SP), H_LN_2_LO, R4 : R4/R7 = N*LN_2_LO 
54 60FD 0740 1 ADDH (RB)+, R4 ; R4/R7 = N*LN_2_LO + LN_FHI_LO 
50 54 60FD 0744 518 ADDH2 R4, RO : RO/R3 = B 
srl a 
are é 1 ; Compute A = N*LN_2_HI + LN_FHI_HI and HLOG(X) 
54 FFO2 CF 8E 65FD 0748 5 § MULH3 (SP)+, H_LN_2_HI, R4 : R4/R7 = N*LN 2 HI 
54 68 60FD O74F 524 ADDH (RB), R4 ; R4/R7 = A = AeCN 2_HI + LN_FHI_HI 
50 54 60FD 0753 525 ADDH R4, RO ; RO/R3 = A + B = ALOG(X) 
05 0757 § $ RSB 
0758 «5 
0798 528 LN_1_PLUS: 
66 11 075 529 BRB LN_1_PLUS_W 
075A 530 
75a «531 
75A 236 ; 
75A 5335 ; Exponent is negative. N =n and F = f 
G73A 534 ;: 
75A 339 
50 57 A2 075A 536 NEG_EXP:SUBW 7, RO ; RO/R3 = F = f 
58 50 07 9C O75D 537 ROTL #7, RO, RB ; R8 = index into MTH$$AB_ALOG table 
58 FFFFFFOO 8F CA 0761 538 BICL #-256 3 = lo exp bit and Ist 7 fract bits 
56 00000000'GF DE 0768 539 MOVAL G*MTHS$$AB_ALOG_V, R6 =: R6 = Address of RTL vector entry 
56 66 C ref 40 ADDL : R6 = Address of MTH$$AB_ALOG table 
58 6648 90 0772 541 MOVE (R6)CRB8), RB : RB = offset into H_FHI fables 
48 19 077 rf BLSS LN_i_PLUS_W ; Branch to special processing 
77 4 ; for F close to 
77 44 
77 45; 
a $$ ; Compute Z, Z**2, P(Z**2) and 2*P(Z**2) 
77 48° 
7E 57 sore 77 49 CVTWH R7, (SP) ; Push N onto the stack 
58 FB7E CFS FD 77 50 MOVAO MTHSSAB_H_FHICR8), R8  ; RB = Address of FHI 
68 63FD 07 31 SUBH (RB), RO, ~R4 + RG/R7 = F = FHI 
§ 60FD 7 é ADDH ( : RO/R3 = F + FHI 
7¢ 54 67FD O78c 8655 DIVHS RO, R4, -(SP) : (SP) 1 = (F = FHID/(F + FHI) 
E 6f $3FD 791 54 MULH3 = (SP), (SP), RO + RO/R3 = 2**2 
FEO3 CF OA SFD 60796 )«=— 555 POLYH RO HLOGLEN2~1, LOGTAB2 : RO/R3 = P(Z*#2) 
0 BE 64FD 79D 36 MULH2 (SP)+, R > RO/R3 = Z*P(2*#2) 
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Oks - Special HLOG/HLOG10 -SEP-1984 91:38: 68 YaTARTE MTHHLOG.MAR; 1 ° ( 
a ; Compute B = N*LN_2_LO + LN_FHI_LO + 2*P(72*Z) 
6 MULH3 (SP), H_LN_2_LO, R4 : R4/R? = N*LN_2_L0 
6¢ ADDH a(R RBS, RG ; R4/R7 = N®LN~ -t0 + LN_FHI_LO 
63 ADDH , RO : RO/R3 = 8B 
65 ; 
66 ; Compute A = N*LN_2_HI + LN_FHI_HI and HLOG(X) 

68 MULHS (SP)+, H_LN_2_HI, R46 = ;_-RG/R7 = N®LN_2_HI 

§ SUBH -(RB), RZ ; R4/R7 = A = AeON 2_HI + LN_FHI_HI 

0 ADDH2 RG, RO : RO/R3 = A + B = ALOG(X) 

71 RSB 

a‘ 

(> 3 

Me ; Special logic for F close to 1 

1 

77 LN_1_PLUS_W: 

78 SOBH3 #1, RO, -(SP) ; (SP) = W=F - 1 

79 POLYH (SP), #LOGLEN1-1,LOGTAB1; RO/R3 = Q(w) 

80 MULH2 | hd RO 3; RO/R3 = WeO(wW) 

81 CVTWH ° : R4/R5 = N 
Ht MULH3 R4, H_LN_2_L0, =-(SP) ; (SP) = N*LN_2.L0 

5 ADDH2 (SP)+> RO ~ ; RO/R3 = N*LR_2_LO + wea(w) 
Se ADDH (SP)+, RO ; RO/R3 = N*LN_2_LO + LN(F) 

5 MULH2 H_LN sul, R4 ; R4/R5 = N*LN-2-HI 
86 ADDH LR 3 RO/R3 = HLOGTXS 
tf RSS 
3 3; X =< 0.0, signal error 
91 ERROR: PUSHL (SP) ; Return PC from JSB routine 
3 MOVZBL #MTHSK_LOGZERNEG, -(SP) ; Condition value 
9 #15, #T, RO ; RO = result = reserved operand -0.0 
94 : Goes to signal mechanism vector 
95 : (CHFSL MCH_RO/R3) so error handler 
36 ; Can modify the result. 
9 CLROQ Rs : 
98 CALLS #2, G*MTHSSSIGNAL 3 Signal error and use real user's PC 
99 3 ependent of CALL vs J 
a1 RSB 5 ate - RO pestered: Re. CHF SL_MCH_RO/R3 
08 
04 END 
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Symbol table ° 6-SEP-1 198 9 £38: $5 Ent THRTL.SRCJMTHHLOG.MAR; 1 . (7) | 1-0! 
ACMASK = 
ERR 4 
— 1 
HLOG = 
HLOG_BASE_10 = 
HLOG_BASE-2 = 
oc C¢ 9 1 
H_INV_LN2_CONS : 8 1 
iH : | 
H-LOGIO E : 1 
LA_1_PLOS 4 
LN-1~PLUS_W 1 
LOGLEN1 = 00000014 
LOGLEN2 = 
LOGTAB1 60 R 01 
LBSI ABE O5A0 R 01 
AT Lf erereree x 00 
MTHS HT 0 RG 01 
MTHSSJACKET_HND tteeeeee Xx 01 
M1H$$S1GNAL eereeeee§ =X 00 
THSHLOG 0000690 RG 01 
MT 000006A6 RG 01 
MTHSHLOG10_R8 00000604 RG 01 
L 444 9d RG 01 
MTHSHLOG RB OO0006E4 RG 01 
MTHSK_LOGZERNEG tereeene 00 
NEG_EXP 0000075A R 01 
= 00000008 
p$ewmwmonr oer woor sce es + 
; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
s OS a 00000000 <¢ 0.) 00 ¢ QO.) NOPIC USR CON ABS LCL — NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE 00000802 ( 2050.) O01 ¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 


Phase Page faults CPU Time Elapsed Time 
Initialization $9 Q: 9? 00: 8 63 
ponnene processing 110 0: . § 0: -64 
s 1 | :00:01. :00:05.48 
oe table sort :00:00.01 :00:00.01 
Pass 2 126 :00:01.39 :00:06.46 
Symbol table output 4 :00:00.04 :00:00.04 
Psect synopsis output :00:00. :00: 2 
Cross-reference output :00: 
Assembler run totals 38 0:00:16.2 


The working set Limit was 1050 pages. 
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11062 bytes (22 pages) of virtual memory were used te buffer_the intermediate code. 
There were 10 pages of synbol table mg > al ecoted F hold 30 non-local and 0 local symbols. 
source Lines were read in Pass 1, duc ing object feateae in Pass 2. 
page of virtual memory was used to def ne 1 macro. 
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! Macro Library statistics ! 


bowser on ene nen ee nee enw ee eoee+ 


Macro Library name Macros defined 
_$255$DUA28:(SYSLIBIJSTARLET.MLB; 2 Tarte | o- 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHHLOG/OBJ=O0BJ$:MTHHLOG MSRC$:MTHJACKET/UPDATE=(ENH$:MTHJACKET) +MSRC$: 
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